Quark Schwinger-Dyson equation in temporal Euclidean space 
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We present an elementary nonperturbative method to obtain Green's functions (GFs) for time- 
like momenta. We assume there are no singularities in the second and the fourth quadrants of the 
complex plane of space momentum components and perform a 3d analogue of Wick rotation. This 
procedure defines Greens functions in a timelike Euclidean space. As an example we consider the 
quark propagator in QCD. While for weak coupling, this method is obviously equivalent to pertur- 
bation theory, for a realistic QCD coupling a complex part of the quark mass and renormalization 
wave function has been spontaneously generated even below the standard perturbative threshold. 
00 ' Therefore, our method favors a confinement mechanism based on the lack of real poles. 
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CD . I. INTRODUCTION 

ly-^ ■ The observable spectra of hadrons represent clear information about the S-matrix at color singlet channels. The 
processes e^e"*" hadrons, r hadronic decays and so on, are significant sources of experimental information in timelike 
regime of momenta. The microscopic description of such phenomena at low energy QCD processes are intuitively 
understood in terms of elementary QCD quanta, although the quantitative description is almost missing. Also other 
quantum field theoretical nonperturbative problems, e.g. confinement, require computation of Green's function (GF) 
r~'' for timelike momenta. For understanding QCD the knowledge of GFs is a crucial matter. As QCD is commonly 
accepted as the ordinary quantum field theory, the amplitudes should be obtainable from elementary QCD GFs which 
^ ' could already encode the information about observables. 

Considering the quark propagator in momentum space, the absence of real poles is a tempting idea of quark 
confinement. The spontaneous dynamical generation of an imaginary part of the quark mass can lead to the absence 
of real pole. We argue that this is the scenario of confinement by showing the model where the imaginary part of quark 
propagator is induced below the expected perturbative threshold. Recall that the threshold value would be otherwise 
uniquely determined just by the quark pole mass. Using the formalism of Schwinger-Dyson equations, we will exhibit 
QQ a realistic scenario, in which the quark mass function, as well as the quark renormalization function, become complex 
for almost all timclike momenta. As usually, the quark propagator remains real for spacelike momenta, where its 
values correspond to the results performed in the standard Euclidean formalism. 

Euclidean space Lattice theory represents the method which in principle provides the information about GFs from 
QQ the first principles. Minkowski space simulation in QCD is not recently feasible because of oscillating phase factor 
in the generating functional. To make the method feasible, the continuation to the imaginary time axis is required 
^ ' and the problem is solved in unphysical Euclidean space. Afterward the continuation back is necessary. Such a 
. ^ , continuation of lattice data to the timelike momentum axis has been performed only very recently [l| providing thus 
' data for very low momenta only. Avoiding large systematic errors, the continuation can be performed only within 
^ [ imposing of an additional global analytical assumption 0. 

In this paper we will solve the quark gap equation which is an alternative way to achieve the non-perturbative 
solution for QCD Greens functions [1, 0] ■ The quark gap equation is the part of the Schwinger-Dyson equations 
(SDEs) which when solved exactly could provide the fully dressed Greens functions as well. As SDEs are an infinite 
tower of coupled integral equations, they require approximation and/or truncation of the SDEs system. Similarly to 
lattice formulation of quantum field theory, the most studies of SDEs are performed in the Euclidean space. The 
trial functions with given analytical properties had been used to make naive continuation to the timelike regime. The 
results of the paper [5,] point towards an analytical structure of the quark propagator with a dominant singularity 
on the real timelike axis, while the nature of this singularity has not been determined with confidence. In the light 
of very recent numerical study the singularity can be a branch point and not a real or complex conjugated poles 
suggested in ^] . The other studies preformed on various assumption also point towards the absence of real pole or 
they at least challenge that the real pole could be a dominant singularity of the quark propagator [1, 0, 9 ■ 

Perturbation theory (PT) is the only known method where such continuation is well understood and massively 
used in practice. In fact in PT at finite order, the analyticity assumptions are more specified: the propagators are 
analytical functions in the whole complex plane up to a real positive semi axis of g^. When PT is reliable, the particles 
are revealed in the GFs poles and branch cuts in momentum space. In such circumstances the tree level single pole 
propagator is dressed within its form constrained by the PT analyticity described above. In general it leads to the 
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known forms of integral representations and dispersion relations for GFs. This has been used in QCD SDEs formalism 
long time ago Q , assuming that the spectral representations remain valid for the full nonperturbative solution. More 
recently, the method of solution based on such spectral representation has been been checked in practice for the 
number of the toy models f^, [lol . [ll| . providing the correct solution only for rather weak coupling, while it appears to 
be inefficient when the couplings exceed certain critical values. It is to be noted, that the position of the branch points 
are uniquely dictated by the mentioned spectral method, which gives us a little freedom for spontaneous generation 
of complex masses. The above mentioned facts do not disprove the spectral method completely, however the practical 
failure of the method could be understood as a sign of weakness of the analytical assumptions. 

In the next Section we propose new approach based on weaker analytical assumption (compared to the PT or 
spectral technique discussed above) and 3d Wick rotation is introduced to "rotate" originally space components to 
imaginary axis. The method is applied to the quark gap equation in the Section 3. The obtained solutions is presented 
in Section 4. 



II. FROM MINKOWSKI SPACE TO TEMPORAL EUCLIDEAN SPACE Et 

In lattice theory and in most of the Schwinger-Dyson equations approaches in the literature the so called Wick 
rotation is used to avoid calculations in Minkowski space, wherein the Green's functions are singular and the integration 
is problematic, especially numerically. Besides the singularities problem, there is another aspect impeding momentum 
integration presented in SDEs, the hyperbolic angle of Minkowskian "spherical" coordinates: 



ko = k cosh 9i 

kx — k sinh 6i cos 62 

ky = k sinh 9i sin O2 cos ^3 

kz = k sinh 61 sin 62 sin 6*3 (2 ■ 1 ) 



runs from - cx) to cx) and most (principal valued here) integrals in terms of these hyperbolic angles cannot be found in 
a closed form. (To that point, there exist a semiperturbative prospect to work directly in Minkowski space, the first 
iteration of SDE was performed in the paper , however single real pole propagator was necessary input to perform 
some integration analytically) 

The mentioned Euclidean space formulation, however, has some drawbacks, since physics involves Minkowski space 
and not Euclidean one. One needs to "rotate back" the results to obtain them for timelike arguments. This rotation 
is basically an analytical continuation on the boundary of perturbative analyticity domain and therefore can be 
frequently ambiguous, especially if the results are numerically obtained. 

In order to circumvent the difficulties stemming from the Minkowski metric or inverted rotation we propose a 
different procedure: instead of Wick rotating the time variable we rotate the space components. Clearly, this way we 
maintain the singularity structure, which for a free propagators stays on the exterior boundary of complex contour, 
but this is a small price to pay for the fact that angular integral are more tractable and especially for the fact that we 
do not need to rotate the variable twice. The method we will present here is used to obtain results at timelike regime 
of fourmomenta, while the correlation functions for spacelike arguments can be evaluated in the standard fashion. 

We assume there are no singularities in the second and the fourth quadrants of complex planes of the complex 
variables kx,ky,kz. Giving the Lorentz invariance, the singularities in the kernels can be functions only of p^, this 
assumption is in agreement with the one used in standard Wick rotation, wherein there are no assumed singularities 
in the first and the third quadrants of the complex fco plane. This happen for instance when the obtained imaginary 
part of the square of the mass function is negative, excluding thus any singularities from the I. and III. quadrants 
and the imaginary ki axis as well. The afore-mentioned Wick rotation is sketched in Fig. [TJ Cauchy theorem gives 
the following prescription for momentum: 

kx,y,z — > ^^1,2,3 , 

d^fc ^ J d'^kEj. , (2.2) 

which in the case of original 3 + 1 is identical with standard Euclidean E "spacelike" one. Note only that, the 
additional i appears when the original Minkowski space is of odd-dimensionality. 
For instance the free propagator of scalar particle then looks 

V--' (2-3) 
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FIG. 1: The Figure shows the integration contour in complex plane for the first space component of the momentum. The 
perturbative singularities (cross) lies in exterior vicinity of the contour. 



with positive square 

=pI+vI+pI+pI (2.4) 

If necessary, one can use the Euclidean definition of Dirac gamma matrices, 7o 74 ; 7 ^IEt ^^.'^ redefined 
gamma matrices satisfy {7^^ , ^'^^ } = 2(5'^" . 

Since the fixed square Minkowski momentum = const hyperboloid with infinite surface is transformed into the 
finite four-dimensional sphere in space, the Cartesian variable are related to the spherical coordinates as usually: 

^ k cos 9 
ki ~ k sin 9 cos f3 
k2 — k sin 9 sin [3 cos 4> 

/c3 = /c sin sin /3 sin (/) . (2.5) 



III. QUARK SDE 



In QCD the quark propagator S is conventionally characterized by two independent scalars, the mass function M 
and renormalization wave function Z such that 

s^) = ^tI^tt; - [^^(^^) - + ' (^-^^ 

or equivalently the functions A^B (where simply M = B/A , A — 1/Z) are used when suitable, noting for the bare 
fermion propagator Sq we have A — 1 and B — wlq. If the function M{p) preserves the real pole in the full propagator 
S{p) then the insertion of Feynman is defines the way of loop momentum integration, otherwise it can omitted. For 
shorthand notation, we will also express the quark propagator in terms of the Dirac vector (Sy) and the Dirac scalar 
{Ss) parts of the propagator: 

S{p)^^Sy{p')^S,{p'). (3.2) 

The gap equation for the inverse of S is 

S-\p) - 5„-i(p)-E(p), 

S(p) = iCaq' J -^T^{q,p)G'^P{p~q)S{qhp, (3.3) 
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where T is the quark-gluon vertex, Ca = TaTa = 4/3 for SU{3) group and G"^ is the gluon propagator, which in 
covariant gauges reads 



where G at tree level reads 



+ 



G(fc')-e 



Jk^ 



(3.4) 



G(fc^) = /c" 



(3.5) 



The gauge parameter dependent term in Eq. (|3.4p remains undressed unless the gauge symmetry is broken, which 
we assume is the case of QCD. 

As QCD is a non-Abelian gauge theory, the GFs are essentially gauge dependent, however in any gauge, the various 
Greens functions are related through the complicated Slavnov- Taylor identities. These constraints are especially 
simplified in the Background Field Gauge (BFG), wherein they simplify to the sort of Ward-Takahashi identities 
01 01 [Hi In this case, the symmetry of the system has not been changed b y g auge fixing procedure. When 
solving Schwinger-Dyson equation the BFG can be usefully further exploited [l^, Il9l . [2g | , even without the exact 
knowledge of the missing and unknown vertices. 

Especially, in BFG the quark-antiquark-gluon- vertex function satisfies QED like Ward identity (WTI): 



fc"r„(p,o = 5-i(p)- 5-1(0 



(3.6) 



where k — p~ I is the gluon fourmomentum. We use the advantage of BFG and we will consider quark propagator in 
BFG in this paper. 



IV. METRIC TENSOR TRUNCATION OF QUARK SDE 



In this section we transform the Minkowski quark SDE into the two-dimensional equation in Euclidean temporal 
space. To proceed this we first specify the approximation of the SDE. In order to go step further beyond the simplest 
ladder approximation we will use the exact WTI of BFG to treat product A: • F in the kernel of SDE (|3.3p . This allows 
to entirely evaluate the contribution wich stem from this term. 

The remaining what needs to be specified is the product of the full vertex F with the metric tensor part g of the 
gluon propagator . As an introductory approximation made in the temporal Euclidean space, we simply take for the 
product g^u^^ ~ lu- The approximation is improvable by making a loop expansion with dressed internal propagators 
(i.e. skeleton expansion). This is a future program which, in addition, will check the reliability of approximation used 
here. 

For convenience we will denote 

S = S]T + Si^ = S<, + Si + Si^ (4.1) 



^^[p)^5B^[p)-5Mp) i> (4.2) 

for i = g,L, L^, where stems from the dressed transverse part of the gluon propagator and L, (Lj) labels the 

selfenergy contribution which follows from the dressed (undressed gauge) longitudinal term in gluon propagator, 
clearly T = g + L in our notation. 

In this notation the appropriate terms explicitly read 

E,(p) = -tZ.g^CA f T^{k,p)g^^''G{q)S{k)Y , (4.3) 

Jk 

Si(p) = iZ.g^CA I r^{k,p)^G{q)S{k)j\ (4.4) 

Jk Q 

= -tZ.g'CA^ lj,{k,p)^S{kh- , (4.5) 
where we have used the shorthand notation J^, for the fourdimensional integral J -A^- 
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Using the WTI we can get for E^^ 



-tZig^CA^ 



(fc2)2 



(4.6) 



The first term (|4.6p is zero since it is odd in the variable k. 

Performing 3d Wick rotation and integrating over the Euchdean angles we get the following contribution to the 
renormalization function; 



SAlAx) 



(47r)2 



Bix) 



JQ 2; 



+ A{x) / dySy{y) 



where x — p\ and y ^ q%. 

For the contribution to the function B we can obtain 



5BlAx) 



2(47r)2 



A{x) / dy!^Ss{y) + B{x) / dyS^iy) 



(4.7) 



(4.8) 



In order to calculate Et the transverse part of the full gluon propagator needs to be specified. At low the 
BFM gluon propagator is unknown function of momenta and gauge parameter, the only known is the undressed 
longitudinal part. To that point we will consider Landau gauge, assuming that various recent studies SDEs and 
lattice calculations performed in this gauge, offer already reasonable estimate. The most ambiguous is the deep 
infrared behaviour q'^Aq^j-,, depending on the details, most of the recent studies shows up that tree level = Q pole 
singularity is softened 'o\ l/{q'^)'^, a < 1, with possible infrared finite solution Tt], T8|. 

In the present paper we assume q'^ is a branch point of gluon propagator, which does not have a purely real pole in 
its transverse part. More specifically, we will assume that the product of the coupling with gluon propagator can be 
expressed through the following integral representation: 



fGiq^AQco)- rd.(^4^^^. (4.9) 

47r ^ Jq q^ - V + i£ 



Thus, contrary to studied quark propagator, the standard analyticity for gluon propagator is still assumed. As already 
mentioned, such a representation has been already used in SDK context [91], however that the exact gluon propagator 
has not the assumed analytical properties is quite possible which would complicate our analysis in this case. 

To do our best we will use the reasonable model of the gluon propagator at all scales. Below, we discuss several 
basic requirements which should be satisfied. 

Firstly, the prescription (|4.9p will respect asymptotic freedom thus for sufficiently large iy^q leading power 
behaviour must be softened by standard perturbative log corrections such that 



' G{q\ Aqco) - ^ , (4.10) 



47r^^'^''-'^^-^-g2Zo5(<z2/A|^^) + . 



where the dots represents higher order scheme dependent contribution. 

It will have no unphysical singularity (known from naive use of perturbative theory at strong coupling). At last 
but not at least, the pg in the gluon propagator may involve confinement. The last two requirements listed above are 
automatically satisfied for any regular function pg. To comply with this we will not assume that pg includes Dirac 
delta as it would be when free particle mode is expected. 

To satisfy all the requirements simultaneously the propagator function G can be constructed by considering the 
function 



where the function a{x) is calculated through 

Pa{x) 



Att/13 



7r2-V {x/Al^^) 



a{x) = P. I dv^^, (4.12) 
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where symbol P. stands for Cauchy principal value integration and j3 in (|4.12p represents the beta function coefRcient, 
for which we take 47r//3 = 1 (recall, 47r//3 — 1.396 for three active quarks in perturbative QCD). 

Recall also, the auxiliary functions pg, a(x) correspond to the imaginary and real parts of the analyticized 1-loop 
effective charge olqcd{—x) constructed in (2ll . [23 , [23| , however the original meaning of a{x) is lost here. In our 
approach it is the gluon propagator, and not the running charge, which satisfies dispersion relation (|4.9p . The full 
expression for a can be found in the original paper. 

Substituting IR (|4.9p into E we can write for 5Ag 



4p2 
-iAttCa 



2p ■ k 

k 



dv- 



where q = p — k. Performing the 3d Wick rotation and integrating over the Euclidean angles we get 

-Ca 



SAgip' 



•7r2 ; 

Jo 



dyy^/y/xSy{y) j dvpg{v)l2{x,y,v) , 



(4.13) 



(4.14) 



where the function I2 is defined below by (j4.17p . Similarly we can easily derive the contribution from g to the function 
B 



dyySs(y) / dvPg{v)I{x,y,v) . 



(4.15) 



The functions /, I2 in (I4.15|) and (|4.14p are the complex non-holomorphic functions defined through the angular 
integral in the following way 



I{x,y,v) = 
l2(x,y,iy) = 



27r 



d0- 



sm 



x + y — — 2^/xycos, 6 + ie 



dO sin OcosO 



Q X + y — — 2y/xy cos9 + ie 



(4.16) 
(4.17) 



Both integrals above can be evaluated in a closed form and we list the results in the Appendix A. 

For L contribution we first use the WTI (|3.6p and the integral representation (|4.9p , then the appropriate contribution 
can be written like 



Jk 9 Jo 



(4.18) 



Making the appropriate trace projections, performing the 3d Wick rotation and after some trivial manipulations we 
get 



5Al{x) = ^ / dvpgiv) 



5Bl{x) 



dvpg{v) 



dyy 
dyy 



~Sy{y)A{x)x 
-S4y)B{x) 



ddq.k sm"^ 9 
q^iq^ ~ V + ie) 
dOa ■ k sin^ e 



SMB{x) 



dOq.psin^e 



where 



q^iq^ - V + ie) 
(27r)3 ■ 



Ss{y)A{x) 



■ p sm" 



q^{q^ ~ v + ie) 
2 



q^iq^ ~ ^ + ie) 



(4.19) 



(4.20) 



In the above formula we do not state explicitly the fact that the all scalar products are in Et space. The scalar 
products k ■ q = k"^ — k ■ p and q ■ p = k ■ p — p^ in the numerators lead finally to the result that can be expressed by 
the integrals / and I2 ■ Explicitly we get: 



5Al{x) 



dyyS^{y)A{x) 
dyyB{x)Ss{y) 



^ dvpg{v) 

V 

dvpg(v) 



\y'i-{x, y, ^) - Vxyhi^x, y, ly) - yl{x, y, 0) + ^/xyhix, y, 0)] 



I{x,y,v) - ^/yjxl2{x,y,v) - l{x,y,0) + y/^y/xl2{x,y,0) (4.21) 
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Similarly the function Bl can be written in the following form: 



6Bl{x) = Cp dyy{Sy{y)B{x)+A{x)Ssiy)) / dvpg{v) 



2 

, \yy{S,{y)B{x)~A{x)Ss{y)) H d,,pg{i.){y - x) ^^'''^'''^ ^ ^^'^'^'^'^ . (4.22) 
Jo Jo 

V. SOLUTION OF SDE IN Et 

Assuming regularity of functions Sg^ Sy on the real axis the quark SDE is transformed into two coupled complex 
integral equations which are free of non-integrable singularities and so they are prepared for suitable numerical 
treatment. Beside, assuming a perturbative asymptotic ultraviolet solution, the SDE requires renormalization. For 
this purpose we use the momentum subtraction renormalization scheme, so the SDE for unrenormalized functions 
A, B which formally reads 

B = mo + ^5B,; A = l + ^M,;i = r,Lj; (5.1) 

i i 

are rewritten into the SDE for renormalized ones. The renormalization constant Zi is absorbed defining thus the 
renormalized propagator, however here we are working in Et space and a certain care is needed. First, we avoid 
the mixture of different computational approaches by choosing a timelike renormalization scale. Further, we keep 
the renormalization constant real, thus only the real parts of the functions SA, SB can be subtracted. Hence the 
renormalization is performed as the follows: 

SAR{p,fi) = ReSA{p) - ReSA{fi) + iImSA{p) , 

SBr{p,p) = ReSB{p) - Re6B{p) + iIm6B{p) , (5.2) 
which leaves us with the renormalized SDE 

Ar{p,h) = 1 + j dy{[ReKA{x,y) - ReKAifi,y)]+iImKAix,y)) , 

Br{p,p) = m{fi) + J dy{[ReKBix,y) - ReKBin,y)]+iImKBix,y)) , (5.3) 
where, for clarity we have explicitly indicate 

j dyKA{x,y) = Y.^^-^ (5-4) 

in order to show how the subtraction procedure works for the integral kernels. The same is performed for similarly 
for the kernel Kb- 

As in the case of perturbation theory, the imaginary part is expected to be finite and untouched by renormalization. 
Clearly, this procedure maintains the hermicity of the Lagrangian. 

For very low momenta the quark masses should approximately correspond to the known values of various constituent 
quark models, where M(0) ~ ^qcd for up and down quarks. Assuming that the real part of the mass function 
is continuous when crossing zero, this value is actually available from Euclidean (spacelike) SDE studies: a typical 
estimate of the infrared mass lays in the range 250 — 600MeT^, while the renormalized mass at few GeV mu,d{'2,GeV) — 
2 — SMeV is the standard input. Here, working in Et space instead of large, we rather choose low renormalization 
scale /X, concretely 

M = Aqcd/4 (5.5) 

adjusting the renormalized function is B{p) = Aqcz? and A{p) = 1. 

In practice the integrals are replaced by the discrete sums on suitable grid. Setting large upper bound c^^Kqcd 
and taking the large number TV = 300 — 1000 of Gaussian mesh points shows up reasonable stability of the numeric. 
The functions A and B are separated to their real and imaginary parts and we solve resulting four coupled integral 
equations simultaneously by the method of iterations. Comparing to the Euclidean spacelike case, the resulting kernel 
of Euclidean timelike SDE (j5.3p is not a completely smooth function a more careful analysis is required. To speed up 
numeric significantly we first integrate over the integral variable v before running the iterations. 
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P/^QCD 

FIG. 2: The quark propagator function A, B. The upper (down) curves represent the real (imaginary) parts. The dimensionfuU 
quantities are rescaled by QCD scale Aqcd- 

Recently we have obtained the results for Landau gauge ^ = 0, where we have achieved a good stability of our 
numerical solution. In Fig. 2 we present the resulting functions obtained for 600 points and e^^Agcu cutoff, enlarging 
cutoff or decreasing the number of points makes the infrared behaviour more chaotic (leaving the smooth average 
approximately constant). 

In PT the propagator is purely real under the threshold scale. Here, this is the main result of our presented study, 
the imaginary parts of the functions B and A are generated below the expected perturbative threshold. The resulting 
mass function becomes complex and a real pole is not present on the real axis of square of momenta. 

More interestingly, we plot the absolute values of the mass function and the inverse of renormalization function in 
Fig. 2. The function \M\ shows up the maximum at 2.3Aqcd where it also cuts the linear function of p. The phase 
(j)M of the mass function extracted from M — |M|e'^'^*'^ where we got cj)M — —25° a,t q — 2. 3 Aqcd- The mass phase 
is a slowly varying function in the full momentum regime and it monotonously goes to small negative value in UV. 

The function A — \A\e^'^'^^ is predominantly real, slowly varying, affecting quantitative behaviour of the function 
M far from the renormalization point. The results for its absolute value and phase are added to the Fig. 2 and Fig. 
3 respectively. Keeping the low scale renormalization point, it has a minimum at few Aqcd and it logarithmically 
increases in UV (the same is also true also for spacelike regime). 

In finite temperature and density QCD it is sometimes suggested, that confinement/deconfinenient phenomena goes 
hand by hand with chiral symmetry breaking/chiral symmetric phases. In our formalism, although there is no space 
for temperature definition, but the description of chiral symmetry breaking could be a QCD must. Nowadays, the lack 
of a reasonably precise description of chiral symmetry breaking is a basic weakness of our presented Et formalism. 
Actually, within our setting, taking the Lagrangian mass to zero (also avoiding forbidden mass subtraction) we got the 
zero dynamical mass everywhere. A bit vaguely pronounced: the kernel of the SDE is not strong enough to produce 
this nonperturbative effect. Without going into technical details, the phenomena of dynamical mass generation could 
be available by further modeling of SDE kernel, (e.g. most naively, by further enhancing of the gluon propagator in 
the infrared). However, as we have found, the price we would pay is an unpleasant (and sometimes drastic) loss of 
numerical stability. 

VI. SUMMARY AND CONCLUSIONS 

We have presented a first analysis of the quark gap equation in the temporal Euclidean space. Given fact that 
3d Wick rotated kernel is non-analytical function at timelike axis of momenta, we do not have at hand the powerful 
method as in the case of standard (spacelike) Euclidean formalism. Nevertheless, at this level the method really works 
and allows us to solve quark gap equation with a good accuracy. 

We obtain the solution with spontaneous infrared complexification of the quark mass function, as opposed to 
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FIG. 3: The absolute value of quark mass function M and the inverse of renormalization function- the function [|j4[| is displayed. 
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FIG. 4; The phase of quark mass function M. The same for the function A. 

the perturbation theory, the quark mass function becomes complex from the beginning of the momentum axis. In 
Landau gauge, B is the main source of the absorptive part of the quark propagator in the infrared region, while the 
renormalization function appears to be marginal for the confinement due the its small generated imaginary part. The 
absolute value of the complex mass function is enhanced at few Aqcd , with the nonzero quark mass function phase 
(j)M — 25° responsible for the absence of the quark propagator pole. 

The method provides not only a qualitative but even a quantitative description of propagator of confined quarks. 
Following the fact that our kernel approximation is only too weak in order to produce correct chiral symmetry 
breaking, we can expect that the observed complcxification phenomena will persist for more realistic kernels of the 
quark SDE. Using the advantage of BFG, the contribution from the longitudinal gluons to quarks selfenergy has been 
already fully taken into account. The product of metric tensor with the improved quark-gluon vertex could provide 
the known slope of the mass function (already known form spacelike studies). To justify our estimate explicitly, an 
improved study of the quark propagator with a stable numeric is required. 
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APPENDIX A: INTEGRAL / 

Consider the integral /: 

I{x,y,v)= de -— (Al) 

Jo a — hcosO + ie 

where a = x + y — v in a, real number and h = 2y/xy is a positive real number. Since v is positive the integrand has a 
singularity in the integration range, so we keep the ie prescription of the propagator. 
Making the standard substitution t = tan | we arrive at the following formula: 

where 

h — a 



c 



b + a 

Performing the principal value integration one can arrive to the following result 



-I{x,y,v)^^ + ^-^-i2^6{c) (A3) 



with an integrable singularity in ^/c, (-v/c) ^ ■ 
The integral 



r/ N sin^^COS^ 

l2ix,y,i^)= de — - A4 

Jo a— ocos^ + ze 



can be evaluated in a similar fashion. 

The result is irregular at c = and regular for positive or negative c. For c > it reads 



l2{x,y,v) 



a + b 



8c 
(l + c)3 



+ 



(A5) 



For c < we can get 



TT 1 - V ^ 

l2{x,y,i') = ————=- , (A6) 
a + b [1 + V-c) 

which is again a finite function. 

The special cases v ^ simplify, the functions l{x,y,0), l2{x,y,0) can be obtained by considering the appropriate 
limits. 
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